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Abstract 



Revisiting canonical integration of the classical pendulum around its 
unstable equilibrium, normal hyperbolic canonical coordinates are con- 
structed. 
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1 Pendulum near the separatrix 

The theory of Jacobian elliptic functions, for reference see [1], yields a com- 
plete calculation for the motion of a pendulum as a function of time. This is 
revisited here, to exhibit a few interesting properties of the elliptic integrals. 
Write the pendulum energy, with inertia moment / and gravity constant 
(rather than the usual g), in the canonical coordinates {B,(3) or as 



where the origin in /? is set at the unstable equilibrium: the definition implies 
that g has dimension of inverse time and the Lyapunov exponents of the 





det 

unstable equilibrium are zizg. B = 1(3 and /? are canonical coordinates for 
the motions. 
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It is well known that near the unstable equilibrium of the pendulum 

5 = 0, /3 = it is possible to define a canonical transformation, mapping 
the origin into itself, introducing new local cordinates (p, q) such that 



B = R,{p,q), /3 = Sc{p,q) (1.2) 

with R,S holomorphic in a polidisk \p\, \q\ < k with k > 0, and in terms 
of which the motion near B = (3 = is described by a Hamiltonian G 
depending on the product pq only, of the form U{p ■ q) = H(B,I3) with 

The purpose of this note is to derive a proof of the latter statement 
via the theory of elliptic functions: this is not the simplest approach if one 
is just interested to know the existence of normal hyperbolic coordinates: 
existence of R, S could be easily established without deriving their "explicit" 
expressions for p and q in terms of elliptic functions. Here we also correct a 
few errors in the earlier attempt made in [6, Appendix 9]. 

The natural correspondence between the hyperbolic fixed point of the 
pendulum and its elliptic fixed point is briefly reported in Appendix B and 
leads to the construction of the normal canonical coordinates for the small 
oscillations, hence to the action-angle variables. 



2 Solution in terms of elliptic integrals 

Motions near the unstable equilibrium have a quite different nature depend- 
ing on the sign of the total energy H{B,P) = U: the ones with U < 
arc "oscillations" (their motions do not encompass the full perimeter of the 
circle) while the ones with U > are "librations" . Therefore it will not be 
possible to find global action-angle coordinates: motions near the separatrix 
(which with our conventions has U = 0) require other coordinates to be 
expressed in a simple way. 

Introduce the variables that appear in the theory of Jacobi's elliptic 
functions 



k' =Vr^, h' = — h = y/l-h'^ 

V 1 + K"' 

U =2g I—., u = t\— 90 = 9: 



(2.1) 



]j 21 ^2/i'K(/i) 

where, K(A;) (1 — k"^ sin^ a) ^o(a. Hence the separatrix has k = +oo 

and U = 0; and the data above the separatrix correspond to ?7 > (or 
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A; > 0). Note that 90 (0) = 9 because K(0) = ^ (as C/ = corresponds to 
k = 00 and h' = 1, h = 0); the following formulae become singular asU ^ 0, 
but the singularity is only apparent and it will disappear from all relevant 
formulae derived or used in the following. 

Other important quantities in the elliptic functions theory are, see the 
references [1, (8.198.1)],[1, (8. 198.2), (8. 146)], 



21 + Vh/ l + 2E^=ie(^)^"' 

where ({k') denotes here what in [1] would be q{k') ([1, (8.146.1), (8.194.2)]). 

In terms of the above conventions wc have, directly from the definitions 
of am, cn, sn, dn (Jacobi's elliptic functions, [1, (8.14)]), and from the 
equations of motion: 

P{t) = 2 am {u,ik), u=*-^, B{t) = 10 dn{u,ik) (2.3) 

([1, (8. 143), (8. 141)]). So that the action B is given as a function of time 

2Iq cn(-i^,h) , , 

([1, (8. 153. 9), (8. 153. 3)]) assuming that initial data are assigned with /? = 0. 

The t dependence of B{t) is naturally expressed via the argument ^ = 
khK.{h) ' second of Eq.(2.4) is used, since kh = h', see Eq.(2.1). This 

explains the role that 

^ l.def _Tr g TT g , . 

^°^^^=^°^2fc^=27I^ ('-^^ 

Eq.(2.5), will play in the following analysis an important role, and it admits 
a rather simple product expansion, [1, (8. 197.1), (8. 197.4)], 

9oi-')=9U{T^) (2-6) 

n=l 

and its logarithmic derivative is 4^^^^ i-x'^" ^° ^^^^ ^''^^^S9o{x') is 
log 04{z, x') ^, where di{z, x'), [2, p.463,489], is a Jacobi's theta func- 
tion. 

It is also convenient to remark that in a motion with energy U it will be 
H{Bit),P{t)) = U = 2g^I^ (2.7) 
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3 Power series representation 



From the theory of ehiptic functions the evolution B(t), /3{t) with any initial 
data above the separatrix {i.e. with /3(0) = and B(0) = II3{0) correspond- 
ing to a given value of h, with U > 0), can be expressed as B{t) = R{-f,S) 
and P{t) = 5(7,(5) with 7 = e^°*,5 = e~^°* and, taking into account [1, 
(8.146.11)], 

i?(7,<5) = -45o/E ^" ^ ^ 1 4^.-1 ^ (3.1) 

n=l '- 

with ^ = ^(/i). Definitions in Eq. (2.1), (2. 2) yield 

go(0 = =5(1 + 7^^ + ...) (3.2) 

which is analytic in by [1, (8.113.1)] near h = 0. 

Eq.(2.2) implies that ^ = A + O(A^) is analytic in A near A = so that 
/i^ = 16A + . . . = 16^ + . . .. Therefore 

go = {l + A^ + 12e + ...)9 (3.3) 

is analytic in ^ near ^ = 0. 

The evolution of if is then a consequence of Eq.(3.1) which leads to an 
expression for S by the remark that R = gQl(jd^ — Sds)S (just expressing 
that B is I times the derivative of /?): namely 

5(T.^)--4s: /,^,L^ (3.4) 

n=l ' 

and, after developing in powers of ^ the denominators and resumming, 

00 

5 =4 E ( arctan(e" S^/^) - arctan(^™ 7V?)) 

^ r- (3.5) 

The first formula reminds of one found by Jacobi which he commented by 
saying that " inter formulas elegantissimas censeri debet", [2, p. 509] (i.e. " 
it should be counted among the m.ost elegant formulae ". 

Note that 50 depends only on ^, see Eq.(3.2), which would be surprising 
if the mechanical interpretation was not taken into account. The Eq.(3.5) 
exhibits the convergence of the map {B, P)< — >■ (^, 7), since ^ < 1 in the region 
above the separatrix: in the latter region Eq.(3.5) provides a convergent 
expansion of the solution. 
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4 Hyperbolic Coordinates 



Motions with initial coordinate /3(0) 7^ also admit a rather simple repre- 
sentation. Remark that all pendulum motions with > (hence different 
from the two equilibria) pass at some time through a phase space point 
with P = 0. If /? is their velocity at that moment we can find a quantity ^ 
such that P,P are given by Eq.(3.5) with 7 = 5 = 1 Therefore they can be 
represented, at least as long as t/" > 0, /? > by introducing the dimension- 
less variables q' = 7-\/f , p' = 5^/^ and allowing 6, 7 to be arbitrary. Then 
the motions will be t ^ (p'e^"*, g'e^^"*) showing that the motion can be 
represented by the following two functions, 



S' =A E (arctan((pV)"*g')-arctan((pV)'"p'))> 

m=0 

n -4J go [ i+((pY)m p>y2 + T+((^^yfr^) 



(4.1) 



m=0 

The motions t — ^ (p'e^°*, g'e~f°*) solve the equations of motion if p',q' {i.e. 
7, 6) are positive. But the equations of motion arc analytic, hence the formu- 
lae Eq.(4.1) together with t (p'e^"*, g'e~^"*), with go = goip'q'), give solu- 
tions of the pendulum equations independently of the sign of p', q', provided 
the series converge. The convergence requires \p'q'\ < 1: which represents 
many data, in particular those in the vicinity of the separatrix. 

The coordinates can be called "hyperbolic" being suitable to describe 
motions near the separatrix (where p'q' = 0) . We also see that time evolu- 
tion preserves both volume elements dBd(3 and dp'dq'; which means that the 
Jacobian determinant q^i'^) must be a function constant over the trajecto- 

def 

ries, hence a function D{x') of x' = p'q'. Note that D{x') has dimension of 
an action. 

It is then possible to change coordinates setting p = a{x')p' , q = a{x')q' 
and choose a(x) so that the Jacobian determinant for (B, f5)< — > (p, q) is = 1. 
A brief calculation shows that this is achieved by fixing 



a" 



{x') = ^ r D{y)dy, (4.2) 
x' Jo 

which is possible for x small because, from Eq.(3.5) and (2.2), it is D{0) = 
32Ig > 0. Therefore the variables, which will have the dimension of a, hence 
of a square root of an action, 

p = p'a{x'), q = q'a{x'), (4.3) 
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have Jacobian determinant 1 with respect to {B,/3) and the map {B,f3)< — >■ 
(p, q) is area preserving, hence canonical. The Hamiltonian Eq.(1.2) becomes 
a fmiction U(x) oi x = pq and the derivative of tlie energy with respect to x 
has to be go{x') (because the p,q are canonically conjugated to B,P). Note 
that X has the dimension of an action, while p, q are, dimensionally, square 
roots of action. 

This allows us to find D{x'): by imposing that the equations of mo- 
tion for the (p, q) canonical variables have to be the Hamilton's equations 

with Hamiltonian U{x)'^= U{x') = H{B,(3) it follows that ^^^^ = go{x'), 

i-^- = 9o{x') or ^ = g,{x'){^{x' a{x'n = 9o{x')D{x') by 

the above expression for a{x'). The just obtained relation, together with 
Eq.(2.2), gives 

D{x') = go{,x')-^^U{x') (4.4) 

which is an explicit expression for the Jacobian q\^i'^}^^ = g'l^'P) — d{p'''q') 
(note that the Jacobian between {B,(3) and {p,q) is identically 1 by con- 
struction). Eq.(4.4) is dimensionally correct because x' is dimensionless so 
that U{x') has the correct dimension {i.e. energy). 

The function U{x') is in Eq.(2.7) where P = by Eq.(2.1), is related 
to x' = ^{h) by Eq.(2.2), so that [1, (8.197.3), (8.197.4)], 



Uix') = 2,^7 J, = 2,^7^ = 32V.' g ^'''^ 

n=l 

To complete the determination of the canonical hyperbolic coordinates it 
remains to find an expression for D[x'),U{x) in terms of the elliptic functions 
to obtain the canonical variable and the Hamiltonian in closed form (rather 
than as power series as done so far). 



5 Determination of the Jacobian. Remarks 



It is remarkable is that the function defined above, hence such that 
d 

dx' 



D{x') = -Arix' a^{x')), seems to be simply 



a^z) = 8I^^go{z), (5.1) 

in a common holomorphy domain, for both sides, around z = 0. This is 
suggested by the agreement of the first 200 coefficients of the expansion of 
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the two sides in powers of z: however this is not a proof and the relation 
Eq.(5.1) holds because it can be seen to be equivalent to an identity on 
elliptic functions, as discussed in Appendix A below. 

Remarks: (1) The expansion of D{x') in powers of x' can be derived from 
Eq.(4.5),(2.6), while that of a'^{x') is obtained from Eq.(5.1) and, again, 
Eq.(2.6). 

(2) It is perhaps natural to guess that the function a{x')'^ should be closely 
related to go{x'); this is a guide to its determination as it becomes, then, 
natural to look for it among the derivatives of go with respect to x'. By 
dimensional analysis all ^'-derivatives of I go have the same dimension as a^. 
Looking also at the derivatives of go as candidates for a'^ is an idea due to one 
of us (PG). This follows a similar line of thought which led to a conjecture 
on the canonical integrability of the "Calogero lattice" , [3] , whose proof was 
discovered in two subsequent works [4] and [5] . 

The relation Eq.(5.1) is equivalent to a notable identity between elliptic 
functions, as discussed in Appendix A below. 

Other peculiarities are, setting 321 g = l,g = 1, 

(1) The function go{x'),U{x'), hence -^go{x'),D{x'), have Taylor coeffi- 
cients in powers of x' which are all positive integers as it follows from the 
relations Eq.(4.5) and Eq.(2.6), while U{x) — x seems to have alternating 
sign Taylor coefficients: 

U{x) -x = 2x2 - Ax^ + 20x^ - 132x^ + lOOSx^ + . . . (5.2) 

where U{x) is obtained by power series inversion of x = x'a{x')'^ and from 
l({x) = U{x') together with Eq.(4.5). 

(2) The function U{x'), energy of the pendulum expressed as a function of 
x', has also the form 



U{x') = 32Igl\p'UAp') + q%'{qWVAp') + Q'UAQ')]=x'fix') (5-3) 

which, remarkably, has by Eq.(4.5) to depend only upon x', and have the 
form x'f{x') for some /. This is not a priori evident, unless the mechanical 
interpretation is kept in mind, from the expressions found for U, V, namely 



oo /2e oo f2e+i 
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(3) Existence of an analytic canonical map integrating, near the hyperbolic 
point, the system with energy Eq.(l.l) into one with Hamiltonian U{pq) = 
gpq + 0{{pq)^) is well known: it can be established without an explicit 
calculation by perturbation analysis, see [6, Appendix A3], for instance . 



A Proof of Eq.(5.1) 

Calling E(A;) = {I - k"^ sin^ a)^da it is E(/i) = hh'^^^ + h'^K{h), see 
[1, (8.123.2], and E{h)K{h') +'E{h')K{h) - K{h)K{h') = f, see [1, (8.122]; 
the latter "Legendre's relation", [2, p. 520], combined with ^ = yields 
the identity 

This can be used to obtain an expression for — ^f;^: keeping in mind x' = 
^ ^ ^' ^ ^ It IS - -T^{^^—dh' WJ^~dhr> wmcn is trans- 

formed into '^^"Sa^' _ l„„™/ / 1 'mh') _ 1 dK(h) . 

Form Eq.(Al.l) it follows, therefore, 



d / TT log x' d I IT log x' 

dh' ^^"^ ^ 2 h'h^K{h)K{h>) dh ^^"^ ^ ~2 hh'^K{h)K{h') ■ ^ ' 

and the corresponding derivatives with respect to h are obtained by multi- 
plying both sides by — . 

To establish Eq.(5.1) consider the relation, 

4-\hh''^}-hK{h)=0 (A1.3) 
dh ^ dh } 

see [1, (8.124.1]. This implies by simple algebra, and keeping in mind that 
^ = — the following identity 

fcK(fc) = h"±[hh'\ - i AK(ft) + ^K(ft))) (A1.4) 

which is a known linear equation, solved by K(/i). This can be rewritten, 
since A = as 
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dhy ' ' dhh'K{h) 

2h — d h- 



Remarking that |4 = J^^, Eq.(A1.5) implies, multiplying both sides by 

2 



d^^^^ ^hrK{K)'dh\ ^ ' ^dhhfK{h^') 
_ 2tt d /hh'^K{h)K{h') d 1 
A 7rKf/i'WKr/i) \ 



(A1.6) 



/i'K(/i)d/iV 7rK(/i')/K(/i) 'dhh'K{h)' 
and by the first of Eq.(A1.2) multiplied by ^ = — this is, using k"^ = ^ 

/ dh . d 1 ,\ 
AdloKx'^d^^'K(/i)^^ 



7r2 



dh h'K{h) dh Ulogx'^dh h'K{h) ' 

vr^ d / I d 1 \ 
'.V d^h'-Klh)) 



(A1.7) 



and multiplying by 21 g'^ 4^ it follows 



h'K{h)dh\ dx'h'K{h)- 

dh 
dx' 



or 2 ^ 1 _ or^g ^ d ^ , d TVg 1 ^ 

and setting a(x')2 '^=^ ^■^'£7^ h'K{h) = ^-^lE'doix') the last relation is 
^C/(x') = go{x')^{x'a{x')^) so that Eq.(4.4) and Eq.(4.2) imply Eq.(5.1). 

B Pendulum at the stable equilibrium 

Prom the above results it is straightforward to find the canonical transfor- 
mation that converts the pendulum Hamiltonian in its normal form around 

the stable equilibrium point. The Hamiltonian is now given by Eq.(l.l) 
with the substitution: g = igg- It is natural to define kg = ik in order to 
use the same set of equations from the unstable case. The system energy is 



around the equilibrium point. 



then Us = 2gj / kj and large values of kg correspond now to small oscillations 
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Finally, it is convenient to define 



h'ik) h{k) '^'^^'^ 



and one finds: 

2 K{hs) (A2.2) 
x'^(hs) =e-^'^^^'^y^^^') = -x'{h) 

where we have used [1, (8.128)]. 

With these conventions, and going through computations similar to the 
ones performed to study the unstable point, the relations found for the latter 
can be converted into the corresponding ones for the equilibrium point. 
In particular, by choosing p' = y/x^cos{gQ^^t) and q' = y^sin(5Q*^t) the 
transformation given by Eq.(4.1) is now: 

oo 

S's =4i E (-1)™( arctan((y2 + g'2)m + ^^z)) 
m=0 

{p'^ + q>'^)r-{p'-iq>) 

l-((p'2 + g'2)m(y_^g/))2^ 

where the relation R'^ = g^^I{p'dqi — q'dpi)S's holds. And the energy can be 
written (see Eq.(4.5)): 

Us{x's) = 2gll- = mglx's J] (^2-4) 

n=l 1 + a^s 

The transformation [B, /3) — > {p'.q') is not canonical. The canonical vari- 
ables, {p,q), can be found by looking for a function as{x's) (which depends 
on the constant of motion x's) such that (p,q) = {as{x'g)p' , as{x'g)q') and the 
Jacobian of the transformation is one. It is, as in the hyperbolic case, 

aliz) = -16I^^gi'\z), (A2.5) 
Finally, the normal form of the Hamiltonian now reads: 

Us{x) = 32IglWi-^) (A2.6) 
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where 



W{z) = z(\-2z-\z^ - 20z^ - 132z'^ - 1008z^ . . .) 



(A2.7) 



which can be compared with the hyperboUc case expression Eq.(5.2): 
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